Extended Gapless Regions in Disordered d x 2_ y 2-Wave Superconductors 
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A generalization of the Abrikosov-Gorkov equations for non-magnetic impurities in unconventional 
superconductors is proposed, including higher harmonics in the expansion of the momentum de- 
pendent gap function and a momentum dependent impurity scattering potential. This model is 
treated within a self-consistent calculation to obtain the electronic density of states, the optical 
conductivity, and the gap function in a two-dimensional d x 2_ y 2-wave superconductor. It is argued 
that momentum dependent scattering from the impurities may lead to extended gapless regions in 
the gap function centered around the nodes of the pure d x 2_ y 2-wave superconductor. The associ- 
ated enhancement of the residual density of states may be responsible for the rapid decrease of T c 
and the increase of the London penetration depth with hole doping observed in overdoped cuprate 
superconductors. 



The pair-breaking role of impurities in d-wave super- 
conductors in the weak coupling limit is well known. Re- 
cently the evidence that such weak coupling descriptions 
apply in the overdoped region of the high-T c cuprates has 
been growing, [jjj However, the experiments that track 
the behavior of the gap function have come to divergent 
conclusions. This raises the question of the evolution 
of the gap function in a d-wave superconductor as a func- 
tion of the pairing interaction strength and the impurity 
concentration. In this letter we examine this question, al- 
lowing for angular dependent impurity scattering poten- 
tials. Our conclusion is that the evolution of the gap func- 
tion depends on the character of the impurity scattering 
potential. In particular, we find that extended gapless 
regions grow as the critical point is reached when the im- 
purity scattering potential is weighted towards forward 
scattering. In this case the scattering processes near the 
gap nodes, which connect regions of opposite sign pair- 
ing amplitude, are of increasing importance and reduce 
the gap in this region. This contrasts with the regions 
around the gap maxima, where the impurity scattering is 
less effective. The resulting evolution of the gap function 
is sketched in Fig. 1. 

In conventional Abrikosov-Gorkov (AG) theory, the 
interaction of electrons with non-magnetic impurities is 
governed by a momentum independent scattering poten- 
tial. U Here we generalize this approach to the momen- 
tum dependent case, with a scattering potential Uk,k'- 
The contribution to the self-energy of electrons scattering 
from this static potential is given by the self-consistent 
T-matrix equation, 

Tk,k'(iw n ) = U kik /(73 +^U k ,qCT3G q (iWn)Tq,k'(kjn), (1) 

q 

where Gk(iw n ) is the Nambu single particle Green's func- 
tion in a superconductor, 



G k (iw n ) 



iw„<7 + A k <7i + £k<73 



(2) 



and the Pauli matrices &\ (<7q = 1) form a complete 
basis in Nambu space. The renormalization of the 
bare Green's function is given by the self energy 
£k(icj n ) = rTk.k(i^n), where T is the concentration of 
impurities. In the limits |uk,k'| — * (Born limit) and 
|uk.k'|^°° (unitary limit) the T-matrix is simplified 
considerably. Hence, these limiting cases have been 
studied extensively in the literature. ||,|| However, the 
strength of the scattering potential is not known from 
first principles, and only comparison with experiments 
indicates that some impurities in high-T c cuprates, such 
as Zn, are closer to the unitary limit. H 
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FIG. 1. Schematic illustration of the momentum depen- 
dence of the gap function in a d x 2_ y 2 superconductor (a) with- 
out disorder (r = 0) , and (b) in the presence of non- magnetic 
impurities (F / 0). 

For computational simplicity, the Fermi surface is as- 
sumed to have a cylindrical shape as shown in Fig. 
1, neglecting effects of the more complicated material- 
dependent band structure observed in high-T c com- 
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pounds. J7J In this case the wave vector k at the Fermi 
surface can be simply expressed as an azimuthal angle, <f>. 
This simplifying assumption of a large and approximately 
cylindrical Fermi surface is in qualitative agreement with 
photoemission and Hall angle experiments on overdoped 
and even optimally doped cuprate samples. [||,|| 

The equations for the single particle Green's function, 

Gj^ 1 (iw n ) = Gj^ 1 (iw n ) — Ek(iwn), have to be solved si- 
multaneously with the zero-temperature gap equation, 
A k . = 2 J " D da;„£ k V k ,k'A k / v /£2 + A£, to yield both 
the amplitude of the gap function and the components 
of the renormalized Green's function at a given impurity 
concentration T. ||] Here, Vk,k' is the effective pairing po- 
tential in the d-wave channel of the same functional form 
as Ak- To describe the angular evolution of the gap func- 
tion, we allow for a higher harmonic, Ak = Ao cos(20) + 
AiCOs(6(/>), and simultaneously include such higher har- 
monics in Vk,k' j which we parametrize in factorized form, 
V k ,k' = V(cos(20) + a cos(6</>))(cos(2</>') + a cos(6<£')), 
where a is the parameter which determines the gap shape 
in the absence of impurities. 

In the following, we analyze how the renormalization 
due to electron-impurity scattering affects the ampli- 
tude and the in particular the shape of the gap func- 
tion, A^ = Ao cos(20)+ Ai cos(6</>). Let us examine the 
Born limit. For the choice of the impurity scattering po- 
tential, Uk,k' = 03 [uo + uicos(0 — (j)')), the AG equations 
take the form 



w n = 0J a + r( 



(I> n (Uo + U 2 COS 2 ' 



^ + A| 



A^ = - r( * 1 V U ) cos(2<//) 



(3) 
(4) 



Here (...) denotes the angular average over the Fermi sur- 
face, and constants arising from the radial momentum 
integration are absorbed in T. Note that A^ depends on 
ui n (r.h.s. of Eq. (4)), and hence the actual shape of 
the renormalized gap can be defined only by analyzing 
the corresponding angular dependent density of states 
(DOS). 

There is a gap equation for each of the two components 
of A^. However, due to our particular choice of the im- 
purity scattering potential only Aq is renormalized, 



A = 2N„V p dWn/ (cos(2^) +a cos(60))A^ 

+ A 2 



(5) 



while the gap equation for the second harmonic gives 
Ai = aA . Thus the renormalized gap function is of the 
form: A^(io) n ) = A (io; n ) cos(20)+ aA cos(6</>). Since 
only the first component of the gap function is reduced 
by impurity scattering, Eq. (4), it is evident that the 
node region widens with increasing impurity concentra- 
tion (Ai/Ao -►1/4). 
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FIG. 2. (a) Amplitude of the gap function in a disordered 
d x 2_ y 2-wave superconductor as a function of the impurity 
concentration in the Born limit. A pure forward-scattering 
impurity potential is chosen (uo = Ui = 1), and the gap 
anisotropy parameter is set to a — 0.2. The gap amplitude 
is plotted relative to its value in the absence of impurities, 
r = 0, and V=1.0. The impurity concentrations are mea- 
sured with respect to Toe, the critical concentration at V=1.0. 
(b) DOS corresponding to the specific impurity concentration 
r/Toc = 0.15, indicated by the dashed line in (a). No is the 
DOS in the normal state, and A o = A(V = 0, V = 1.0) (c) As 
(a), but for unitary impurity scattering with uo = 1, Ui = 0, 
and a = 0. (d) As (b), but for the parameters used in (c). 

In Fig. 2(a) the gap amplitude Ao obtained from the 
numerical solution of Eqs. (3-5) is shown as a function of 
impurity concentration and of the strength of the pairing 
interaction. Here we have chosen a=0.2, uo = 1, and ui 
= 1, i.e. the case favoring forward scattering. The su- 
perconducting gap can be destroyed both by increasing 
the impurity concentration beyond a critical value, T c , 
or by reducing the pairing potential V below a critical 
V c . While the first case naturally occurs when introduc- 
ing disorder into samples by impurity doping, the latter 
case may be realized by introducing additional holes into 
overdoped cuprate superconductors. fici[| 

From an expansion of Eqs. (3-5) about T c we find that 
for large impurity concentrations the amplitude of the 
gap function vanishes as Ao cx y/l — T/T c . Similarly, if 
one drives the system critical by reducing V/V c , we find 
Ao oc y/V/Vc — 1, in agreement with the full numerical 
solution of these equations. 

From an analytic continuation of Eqs. (3) and (4) onto 
the real frequency axis, the electronic DOS can be eval- 
uated as N(w) = -| Ek Im [Gk(iwn)]|i<i n =<i+i<5- In Fig. 
2(b) the DOS is shown at r/r c = 0.15 (indicated by the 
dashed line in Fig. 2(a) ). As the strength of the pairing 
potential is lowered, the system approaches the critical 
region where the gap is small, and hence the residual 
DOS increases with decreasing V. As will be discussed 
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in the following, this increase in N(0) is enhanced by the 
occurrence of flat regions around the nodes. 

The above results were obtained in the Born limit. In 
the unitary limit, the AG equations do not reduce to sim- 
ple expressions like Eqs. (3-5) if the impurity scattering 
potential is taken to be momentum dependent. Hence, 
instead of attempting to solve the AG equations for the 
general Ukk' used above, we restrict ourselves here to the 
case Uo = 1, ui =0, i.e. only isotropic impurity scat- 
tering. The results are shown in Fig. 2(c) and (d). It is 
obvious that in the limit of strong scattering the gap am- 
plitude is reduced more rapidly by the introduction of im- 
purities than in the Born limit (comparing Fig. 2(a) and 
(c)). When approaching the critical regime by decreasing 
V/V c , a reduction of the gap amplitude (position of the 
peak in N(w)) and an increase in the residual DOS is ob- 
served, similar to the Born limit treated above. However, 
since only isotropic impurity scattering was considered in 
the latter case there is no widening of the node regions in 
the gap function and hence no additional enhancement 
of N(0). While from the experimental side it has not 
been settled what the strength of the effective scattering 
potential, |u klk '|, should be, a comparison of the two ex- 
treme limits suggests that no dramatic differences are to 
be expected within AG calculations. Qualitative dif- 
ferences between these two limits do appear for the single 
impurity problem, where anisotropic impurity resonances 
appear in the unitary limit. fl2]| 

Let us now turn to the qualitative changes in the shape 
of the gap function in the presence of a momentum de- 
pendent impurity scattering potential. In the clean case, 
there are four nodes located at 4> — (2n~l)ir / A (n=l,...,4) 
(Fig. 1(a)). The first of these nodes is shown in Figs. 3(a) 
and (b) for a fixed impurity concentration (r/r c = 0.15) 
and various values of V. In Fig. 3(a) only isotropic im- 
purity scattering was considered (uo = 1, Ui = 0) while 
in Fig. 3(c) the case favoring forward scattering is shown 
(uo = ui = 1). In the latter case the gap function was 
extracted by analyzing the angle resolved DOS: A(</>) was 
defined as the position of the inflection point in N(^>, ui) 
(maximum in dN((f>,u))/du)). In both cases the gap am- 
plitude decreases with decreasing V. (A similar behavior 
occurs when V is kept fixed and T is increased.) In the 
case of forward scattering, however, a flattening of the 
gap function around the nodes is observed. To quantify 
this behavior, we plot the extension, d^, of the "flat" part 
of A((/>) in the insets. With d^ we denote those segments 
of A(</>) where the magnitude of A(</>) has fallen below 
half of its maximum value, d^ = tt/2 — 20(A max /2). ]T^ ] 
Note that as long as only a finite number of harmonics 
in the expansion of A(0) is considered, only extended 
saddle points can occur. However, it is obvious from 
the above discussion that there is a clear tendency to- 
wards gapless ("normal") areas in the presence of mo- 
mentum dependence in the impurity scattering potential 
coexisting with gapped regions (as depicted schemati- 
cally in Fig. 1(b)). This prediction may be verified in 
overdoped cuprate superconductors by tunneling experi- 



ments, and by angle resolved photoemission spectroscopy 
(ARPES) if an adequate energy resolution (of order 1 
meV) can be achieved. B 
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FIG. 3. (a) Angular dependence of the magnitude of the 
gap function in the Born limit at r/To c = 0.15, uo = 1, 
Ui = 0, and a — 0.2. The inset shows the length of the 
segment of the gap function where |A(0)/A(O)| < 0.5. (b) 
Optical conductivity corresponding to the parameters chosen 
in (a). (Too is the normal state value of the optical conduc- 
tivity, (c) As (a), but with ui = 1. (d) As (b), but with 
ui = 1. 

In the critical regime, where the amplitude of the 
gap function is small, Eqs. (3-5) can be linearized. 
While the renormalization of Ao remains frequency de- 
pendent in this limit, a functional form of the renor- 
malized gap function can be extracted from the angu- 
lar dependent residual DOS, N(0, 0). We find that A = 
Ao[l + Vl\ /(4u§ + 2uf)] _1 , while Ai remains unrenormal- 
ized. For the above choice of parameters (uo = Ui = 1, 
a = 0.2), we then find Ai/A » 0.233. 

A complementary experimental probe that is sensitive 
to extended flat regions in the gap function is microwave 
spectroscopy, which measures the real part of the c-axis 
infrared optical conductivity, 

7re 2 f u f d 2 k 
a{u) = —J^ dui'J — ptr[A k (a,- W ')A k (-u/)], ( 6 ) 

where A k (w) = -ilm[G k (icI> n )]| i£jn=ti+ i,5, and an angular 
average has been taken over the in-plane directions. |l4|] 
In Fig. 3(b) and (d), the optical conductivities are shown 
for the two cases of isotropic and forward impurity scat- 
tering. In the clean case, cr(ui) shows a "knee" feature 
at 2A 00 (A 00 = A(r = 0,V = 1.0)). In the presence of 
impurities, this energy scale is reduced when approach- 
ing the critical regime by decreasing V. In analogy to the 
discussion of the DOS, the increase of <r(0) is enhanced 
considerably in the case of forward scattering. 
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In Fig. 4(a) and (c) the residual DOS is shown for 
the Born limit at uo = 1, ui = 1, and a — 0.2, and 
for the unitary limit at uo = 1, Ui = 0, and a = 0. 
By expanding about the clean limit, it is seen that in 
the Born limit, the zero-frequency DOS vanishes expo- 
nentially slowly, N(0) cx exp(— const/r), for small impu- 
rity concentrations. On the other hand, in the unitary 
limit an infinitesimal amount of impurities is sufficient 
to yield a finite residual DOS, N(0) cx %/T- Hence an 
enhancement of the residual DOS can be induced both 
by increasing the strength of the scattering potential and 
by introducing angular dependence in it. However, only 
the latter case leads to extended flat regions in the gap 
function. 

As a consequence of the finite residual DOS, the low- 
temperature London penetration depth has the form 
A(r, T) = A(0,0)+ N(r,0)/2+ 0(T 2 ). H As the effec- 
tive electron-electron interaction strength is decreased, 
i.e. when introducing additional holes into optimally 
doped or overdoped cuprate superconductors, A grows 
rapidly. This behavior, illustrated in Figs. 4(b) and (d), 
is in agreement with recent measurements by Locquet et 
al. Il6| 
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FIG. 4. (a) Residual DOS as a function of impurity con- 
centration in the Born limit. Here uo = 1, Ui = 1, and 
a = 0.2. (b) Low-temperature London penetration depth as 
a function of the e~ — e~ pairing potential. The parameters 
are chosen as in (a), (c) As in (a), but in the unitary limit, 
and with uo = 1, Ui = 0, and a = 0. (d) Low-temperature 
London penetration depth for the parameters chosen in (c). 

In our results the maximum value of the gap is reduced 
consistently as either T is increased or V is reduced. This 
agrees with the results of ARPES and tunneling exper- 
iments. It disagrees with the conclusion from studies of 
the optical conductivity. However, there is an increase in 
the gapless region in the case of forward weighted impu- 
rity scattering which goes some way towards explaining 
the optical experiments. Our conclusion is that varia- 



tions in impurity character may be at least partly re- 
sponsible for the divergence of the experimental results, 
although an explanation for the absence of renormaliza- 
tion of the gap maximum is not possible within this weak 
coupling scheme. 
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